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ABSTRACT 

I propose a feedback model to explain the correlation between the supermas- 
sive black hole (SMBH) mass and the host galaxy bulge mass. The feedback is 
based on narrow jets that are launched by the central SMBH, and expel large 
amounts of mass to large distances. The condition is that the jets do not pen- 
etrate through the inflowing gas, such that they can deposit their energy in the 
inner region where the bulge is formed. For that to occur, the SMBH must move 
relative to the inflowing gas, such that the jets continuously encounter fresh gas. 
Taking into account the relative motion of the SMBH and the inflowing gas I de- 
rive a relation between the mass accreted by the SMBH and the mass that is not 
expelled, and is assumed to form the bulge. This relation is not linear, but rather 
the SMBH to bulge mass ratio increases slowly with mass. The same mechanism 
was applied to suppress star formation in cooling flow clusters, making a tighter 
connection between the feedback in galaxy formation and cooling flows. 



1. INTRODUCTION 

The tight correlation between the supermassive black hole (SMBH) mass, Mbh, and the 
velocity dispersion, a, of the hot component of the host galaxy (e.g., Merritt & Ferrarese 
2001; Gebhardt et al. 2000; Giiltekin et al. 2009 ) is well established. A recent paper by 
Giiltekin et al. (2009 and earlier references therein) derives the relation as 

log(M BH /M ) = (8.12 ± 0.08) + (4.24 ± 0.41) log((x/200 km s" 1 ). (1) 

There is still no consensus on the slope, though. Shen et al. (2008), for example, find the 
slope in luminous broad-line AGN to be ~ 3.1 — 3.6 instead of 4.24. 

The correlation between the SMBH mass and the the bulge mass, Mb u i ge (e.g., Kormendy 
& Richstone 1995; Laor 2001), is less tigh. Using the Mbh oc £ L11±0 - 18 correlation as derived 
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by Giiltekin et al. (2009), where L is the bulge stellar luminosity, and substituting L from 
the M bulge oc L 118 relation from Magorrian et al. (1998), one derives M B h oc M^ ulge with 
(3 = 0.94 ± 0.15, which is compatible with a linear relation (3 = 1. However, Laor (2001) 
derives a value of (3 ~ 1.2 — 1.8, depending on the exact samples and relations used. Laor 
then deduced that there is no linear relation between Mbh and Mb u i ge , but rather that the 
ratio Mmi/Mbuige increases from ~ 0.05% in the least luminous bulges (My ~ —18) to 0.5% 
in the bright ellipticals (My ~ -22). 

To further explore this relation, I examine the correlation among the 20 brightest and 
20 faintest galaxies in the sample of 31 galaxies used by Giiltekin et al. (2009) in their figure 
4. I find that for the brightest galaxies the value of 7 in the relation Mbh oc L 1 is larger by 
~ 25% than the value of 7 for the entire sample, while for the 20 faintest galaxies it is lower 
by ~ 14%. Using then the relation Mb u i gc oc L (Magorrian et al. 1998), one finds that 
the ratio MBH/Mb u i ge increases with mass, as claimed by Laor (2001), and obtained in some 
models and simulations (e.g., Shabala & Alexander 2009). 

The relation L — a is not tightly constrained either. Laor (2001), citing Nelson & Whittle 
(1996), takes L ~ a 3 , and using M bu ig e oc L 118 (Magorrian et al. 1998) gets M bu i ge cx~ ct 35 . 
Using this relation in equation ([1]) gives (3 = 1.2. Haring & Rix (2004) find (3 = 1.12 ± 0.06, 
and Mbh /Mbuige = 0.14% ± 0.04% for M bu i ge ~5x 1O 1O M . 

Based on these considerations, and for the purpose of the present paper, I assume the 
deviation from a linear Mbh — Mb u i ge relation to be due to an additional dependance on the 
dispersion velocity 



where O ~ 10~ 3 and \ ~ — 1. The above expression is an average one, as the relations 
between the SMBH mass, the bulge mass, and the dispersion velocity are somewhat different 
for elliptical galaxies, classical bulges and pseudo-bulges (Gadotti & Kauffmann 2009). In 
this first presentation of the idea I ignore these differences. 

In this paper I will try to account for relation (J2]) with a feedback mechanism based 
on jets launched by the SMBH. The feedback mechanism where AGN jets (outflow; wind) 
suppress gas from cooling to low temperatures and from forming stars was discussed for 
both cooling flows in galaxies and clusters of galaxies (e.g., Binney & Tabor 1995; Nulsen 
& Fabian 2000; Reynolds et al. 2002; Omma & Binney 2004; Soker & Pizzolato 2005), and 
in galaxy formation (e.g., Silk & Rees 1998; Fabian 1999; King 2003; Croton et al. 2006; 
Bower et al. 2008; Shabala & Alexander 2009). I note that some of the papers (e.g., Silk 
&; Rees 1998; King 2003) make use of the Eddington luminosity limit; the model proposed 
here make no use of the Eddington luminosity limit. Most models (e.g., Silk & Rees 1998; 
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Fabian 1999) do not consider the geometry explicitly; here the geometry of the narrow jets 
and the motion of their source are key issues. Bower et al. (2008) numerically tried to derive 
the M B h — ^buige correlation. In the present paper I try to present the basic physics that I 
propose leads to this correlation. 

2. THE JET-FEEDBACK MECHANISM 

I do not consider the formation of the seed SMBH, but rather assume its existence. The 
proposed mechanism is based on the following assumptions. 

1. The feedback mechanism, i.e., the one that expels mass from the inner regions to large 
distances and heats it up before stars can be formed, is driven by jets launched by an 
accretion disk around the SMBH. This is supported by the presence of X-ray deficient 
bubbles in clusters of galaxies that show that AGN jets are sufficiently energetic to 
heat the intracluster medium (see review by McNamara & Nulsen 2007). 

2. The properties of jets launched by SMBH have some universal properties, such as the 
fast jet's speed Vf ~ c. 

3. There is a universal efficiency (e.g., Tremaine 2005) of converting accretion energy to 
kinetic energy of the two jets. Considering different works (e.g., Kording et al. 2008; 
Merloni & Heinz 2008) and taking into account the result of Binney et al. (2007), I 
take this ratio to be ~ 0.05. As the jets are launched at velocities of Vf ~ c, I take the 
ratio of mass launched in the two jets to accreted mass to be r\ = Mj/M acc ~ 0.05. 

4. The mass flowing in at early stages, i.e., the mass available for star formation, is very 
large. Namely, the mass that is converted to stars is limited by the feedback mechanism 
and not by the mass available in the SMBH surroundings. This is supported by studies 
of galaxy formation (e.g., Bower et al. 2008). 

5. There is a relative motion between the SMBH and the inflowing mass of i> re i — cr, as the 
relaxation time of the SMBH is longer than the galaxy formation time scale (Tremaine 
2005). 

6. The surrounding mass M s that resides at a typical distance r s and having a density 
p s , (see below) is flowing inward at a velocity of ~ a. Thus, M s ~ A-Kr 2 s ap S) and it is 
resupplied on a time scale of ~ r s /a. This assumption should be better constrained 
by future 3D numerical simulations. 
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With these assumptions in hand, I proceed to describe the interaction of the jets with 
the surrounding inflowing gas. If the jets penetrate through the surrounding gas they will 
be collimated by that gas, and two narrow collimated fast jets will be formed, similar to the 
flow structure in the simulations of Sutherland & Bicknell (2007). By fast it is understood 
that the jet's velocity is not much below its original velocity. If, on the other hand, the 
jets cannot penetrate the surrounding gas they will accelerate the surrounding gas and form 
SMW (slow-massive- wide) outflow (Soker 2008). The conditions for the jets not to penetrate 
the surrounding gas but rather form a SMW outflow are derived below. The derivation 
will not include relativistic effects, although the fast jets might be relativistic. The other 
uncertainties involved are larger than those introduced by this assumption (assuming the 
jets are not highly relativistic). 

Let the jets from the inner disk zone have a mass outflow rate in both directions of 
Mf, a velocity Vf, and let the two jets cover a solid angle of 4n5 (on both sides of the disk 
together). The density of the outflow at radius r is 

M f 

Pf ~ 4^V (3) 

Let the jets encounter the surrounding gas residing within a distance r s and having a typical 
density p s . The head of each jet proceeds at a speed Vh given by the balance of pressures 
on its two sides. Assuming supersonic motion this equality reads p a v\ = Pf(vf — v h ) 2 , which 
can be solved for v h 

\ 1/2 / \ 1/2 

AirSrgVfPs \ I 8M s Vf 




M f I \ Mi 



a 



1225 | MflMl] (A.) rJl) 1/2 ( ° Y 112 . (4) 
1 10 4 J \0.1 J \cJ V200 kms- 1 / v; 

where in the second equality the mass inflow rate M s ~ A7rp s ar 2 (by assumption 6), has 
been substituted. The time required for the jets to cross the surrounding gas and break out 
of it is given by 

/ ■ \ 1/2 

r s r s / 8M s Vf 

r^j I J 



v h Vf y MfO J V 1 k P c / 

where in the last equality the same values as in equation (j3J have been used. 

If there are no changes in the relative geometry of the SMBH and inflowing mass, the 
jets will rapidly penetrate the surrounding gas and expand to large distances. In this case the 
jets will not deposit their energy in the inflowing gas. For an efficient deposition of energy 
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to the inflowing gas, we require that there will be a relative motion between the SMBH and 
the inflowing gas, such that the jets continuously encounter fresh mass. The relevant time 
is the time that the transverse motion of the jet crosses it width r s = Dj/v rc \ ~ Dj/a, as by 
our assumption 5 the relative velocity is v ie \ ~ a. The width of the jet at a distance r s from 
its source is Dj = 2r s sina, where a is the half opening angle of the jet. For a narrow jet 
sin a ~ a ~ (25) 1 ' 2 , and 

2(25) 1 / 2 r s 



a 

The demand for efficient energy deposition, t s < t p , reads then 



(6) 



M s ^ v f 

Z 8- 
M f a 



>8^. (7) 



This result can be understood as follows. The ratio Vf/cr comes from the ratio of the ram 
pressure of the narrow jet to that of the ambient gas which disturbs the jet, and from the 
relative transverse motion of the jet and the abient gas. The number 8 comes from the 
geometry of a narrow jet with a relative transverse velocity to that of the ambient gas. 
Using assumption 3 that Mf = r]M acc , and substituting typical values, this condition reads 



L &cc 



> ^ = 600 (i) (^) ^ ^ 



The accretion rate M acc is the accretion rate onto the SMBH, and the inflow rate of 
the surrounding gas is assumed to form stars in the bulge (if it is not expelled by the jets). 
If the inflow rate is above the value given in equation (jSJ), the deposition of energy by the 
jets is efficient enough to expel the mass back to large distances and heat it (Soker 2008). 
Following assumption 4 above, I take the SMBH to bulge mass ratio to be equal to M acc /M s . 
Equation (JE1) yields 

M BH = 0.0017M bulEe ( -^L-Y 1 ( ( -V (9) 

g V0.05/ \cJ V200 kms- 1 / KJ 

The last equation has the form of equation (j2J), and by using typical properties of AGN 
jets and bulges (77; Vf] a), the numerical value of equation (j2j) is reproduced with x = 1. 
Basically, the equation contains one parameter p m = rjVf. Its physical meaning is that of the 
momentum of the material ejected in the jets per unit accreted mass to the SMBH. For a 
more detailed comparison to observations, future work is needed to incorporate the proposed 
feedback mechanism into numerical simulations of galaxy formation and evolution, e.g., such 
as those described by Di Matteo et al. (2005). 



- 6- 



3. DISCUSSION AND SUMMARY 

I built a simple feedback mechanism where jets launched by the central supermassive 
black hole (SMBH) interact with the surrounding inflowing gas in a young galaxy. By using 
several very plausible assumptions, most of which are well supported by observations, I 
derived a relation between the mass that is converted to stars and the mass accreted by the 
central SMBH (eqs. [8] and [9]). The basic process to suppress star formation requires that 
the jets do not penetrate the surrounding gas; if they do, their energy is deposited at large 
distances, rather than in the gas that potentially can form stars in the bulge. For that, the 
jets should encounter new material before they escape. Namely, the typical time for the jets 
to cross its own width at radius r s should be shorter than the penetration time at radius r s , 
t s < t p . This leads to equation (jSJ). 

If the mass inflowing rate M s is larger than the value given by equation (jSj) the jets 
are very efficient in depositing their energy to the inflowing gas. The jets are shocked, and 
since the radiative cooling time is long (Soker 2008), there is no energy loss and the jets can 
expel huge amounts of mass back to large distances and heat the gas. This suppresses star 
formation. Even the Compton cooling time (e.g., King 2003) is longer than the jet flow time 
T a jvf for r s > 10 pc, and here r s ~ 0.1 — 1 kpc. I therefore take the mass inflow rate into the 
bulge to be as the limiting value in equation (|S]). This leads to equation (JSJ, which has the 
form and numerical value as in equation (j2J). In addition to the assumptions of the model, 
equation (jUJ) contains only one parameter p m = r]Vf, which is the momentum of the material 
ejected in the jets per unit accreted mass to the SMBH. 

I note that the criterion that the jets do not penetrate the surrounding material might 
be relevant also to models that attribute the correlation between the SMBH mass and the 
host galaxy bulge mass to a negative feedback (e.g., Silk & Norman 2009); in a negative 
feedback model the jets induce star formation in the surrounding gas. 

King (2003) compares the cooling time to the slow outflow time. Since the outflow gas 
has a velocity < a in his model, he finds cooling to be very important. King assumes the 
presence of an almost spherical shell, and does not consider its formation by two narrow 
jets as the present model does. Once the condition found here is met, i.e., the jets do not 
penetrate the surrounding gas, a slow expanding shell is formed, and some of the calculations 
of Silk & Rees (1998), Fabian (1999), and King (2003) are applicable here as well. However, 
here I made no use of the Eddington luminosity limit, and a similar mechanism can be 
applied to cooling flows (see below). 

As the system evolves it relax. The relative motion of the SMBH and the inflowing gas 
decreases with time. Therefore, as the system ages, and masses and a increase, condition 
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(JHJ) is more difficult to fulfill, and the accretion rate to the bulge can be larger than given 
by equation (jSj). Namely, the SMBH mass is lower than given by equation (jUJ). This implies 
that the power of a there is \ < 1, as is deduced from observations in equation (j2j). 

Although the proposed feedback mechanism has several assumptions, they are all very 
plausible. The feedback mechanism will have to be studied by 3D numerical simulations, 
as the relative motion of the SMBH and the inflowing gas must be included. There is a 
big advantage to the proposed mechanism, in that the same basic mechanism can work 
to suppress star formation in cooling flows (Soker 2008). Using the same mechanism with 
similar parameters, I developed a mechanism for the formation of slow massive wide (SMW) 
outflows (or SMW jets) in cooling flow clusters (Soker 2008). Such jets can inflate the 
observed X-ray deficient bubbles that are observed in clusters of galaxies (Sternberg et al. 
2007). 

I thank Craig Sarazin for very helpful comments. This research was supported by the 
Asher Fund for Space Research at the Technion, and the Israel Science foundation. 
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